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Abstract. The internal gravity waves of low frequency 
which are emitted at the base of the solar convection zone 
are able to extract angular momentum from the radiative 
interior. We evaluate this transport with some simplify- 
ing assumptions: we ignore the Coriolis force, approxi- 
mate the spectrum of turbulent convection by the Kol- 
mogorov law, and couple this turbulence to the internal 
waves through their pressure fluctuations, following Press 
(1981) and Garcia Lopez & Spruit (1991). The local fre- 
quency of an internal wave varies with depth in a differ- 
entially rotating star, and it can vanish at some location, 
thus leading to enhanced damping (Goldreich & Nichol- 
son 1989). It is this dissipation mechanism only that we 
take into account in the exchange of momentum between 
waves and stellar rotation. The flux of angular momentum 
is then an implicit function of depth, involving the local 
rotation rate and an integral representing the cumulative 
effect of radiative dissipation. We find that the efficiency 
of this transport process is rather high: it operates on a 
timescale of 10 7 years, and is probably responsible for the 
flat rotation profile which has been detected through he- 
lioseismology. 

Key words: Hydrodynamics, turbulence; Sun: interior, 
rotation; Stars: interiors, rotation. 



1. Introduction 

Although it is well known that waves do carry momen- 
tum, this process has received little attention so far in 
stellar physics. The transport of angular momentum by in- 
ternal waves (also called gravity waves) has been studied 
in the context of tidal braking involving massive binary 
stars (Zahn 1975a; Goldreich & Nicholson 1989), but only 
recently has it been invoked as a mechanism which could 
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shape the rotation profile in the Sun (Zahn 1990; Schatz- 
man 1993). In contrast, the importance of the momentum 
transport by such waves has been recognized long ago in 
atmospheric sciences (cf. Bretherton 1969a, b): this phe- 
nomenon is responsible in particular for the so-called clear 
air turbulence. 

The purpose of the present paper is to assess the ef- 
ficiency of angular momentum transport in a solar-type 
star, where such waves are generated by the turbulent 
motions of the convection zone. In this first approach we 
shall make several simplifying assumptions, fully aware 
that the outcome must be considered as a crude approxi- 
mation. The results may be easily transposed to massive 
stars, where these waves are emitted by the convective 
core. 

The reason for examining the role of such waves in 
the Sun is that the other mechanisms which have been 
analyzed so far seem unable to achieve the flat rotation 
profile revealed by helioseismology (cf. Brown et al. 1989). 
Both rotation-induced turbulent diffusion (Endal & Sofia 
1978; Pinsonneault et al. 1989) and wind-driven meridian 
circulation (Zahn 1992) fail to extract sufficient angular 
momentum from the radiative interior (Chaboyer et al. 
1995; Matias & Zahn 1996). Magnetic torquing looks at 
first sight more promising, but the field lines anchored in 
the differentially rotating convection zone would probably 
enforce non uniform rotation below (cf. Charbonneau & 
MacGregor 1993), and this is not observed. 

We begin by recalling the main properties of the in- 
ternal waves (§2), and calculate the flux of angular mo- 
mentum carried by a monochromatic wave (§3). We then 
deduce the energy spectrum of those waves from their 
coupling with the turbulent motions at the base of the 
convection zone (§4) and integrate the angular momen- 
tum flux over the whole spectrum (§5). We finally derive 
an estimate for the efficiency of the angular momentum 
transport by such waves in the Sun (§6) and end by some 
concluding remarks. 



2 



J.-P. Zahn, S. Talon and J. Matias: Angular momentum transport by internal waves 



2. Properties of internal waves 

The properties of internal waves propagating in stellar in- 
teriors have been described by Press (1981) and by Gold- 
reich & Nicholson (1989). Let us recall the main results, 
and adapt them to a differentially rotating star. We use 
the spherical coordinates (r, 9, (f>) suited for this problem; 
e z is the unit vector on the rotation axis, and that in 
the azimuthal direction. We further assume that the angu- 
lar velocity f2 depends only on depth, because differential 
rotation in latitude is severely limited through hydrody- 
namical instabilities (Zahn 1975b, 1992). The velocity field 
with respect to an inertial frame is then 



U(r, 9, <P, t) = Q(r)e 2 xr + u(r, 9, <j>, t), 



(1) 



with u being the velocity associated with the wave. The 
equation of motion reads 



^ + |2£le z x u + e ( prsin9u ■ Voj 



p p 



with 

dt ~ \dt + dcj))' 



(2) 



(3) 



and the usual notations for the gravity g, and for the Eu- 
lerian perturbations of pressure (P') and density (//). We 
have simplified this equation by neglecting the fluctuations 
of the gravitational potential: the Cowling approximation 
is amply justified here, considering the high radial order 
of the waves. Viscosity also has been ignored. We add the 
continuity equation 



(4) 



and the energy equation, which in the adiabatic limit re- 
duces to 
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u r = 0, 



(5) 



Ti being the adiabatic exponent. 

We now proceed with a further simplification, which 
admittedly is much less justified. Neglecting the terms in 
curly brackets in the equation of motion (2), we treat the 
waves as if they were pure gravity waves which are not 
modified by the Coriolis acceleration, but just feel the en- 
trainment by the differential rotation. 

Then the differential system above is very similar to 
that governing the internal waves in a non-rotating star; 
the only difference lies in the derivative with respect to 
time, which is replaced here by (3). The solutions are sep- 
arable in spherical functions and time, as shown for the 
vertical component of the velocity: 



where 

a(r) = cr — mO(r), 



(7) 



with <jq being the frequency in the inertial frame. In the 
low frequency range of the internal waves which are con- 
sidered here, the function ^f(r) = p^r 2 u v obeys the second 
order equation (cf. Press 1981): 



d 2 * 
dr 2 
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(8) 



The buoyancy (or Brunt- Vaisala) frequency N is given by 



N 2 = N 2 + N 2 = |^(V ad - V) + |^V M 
lip rip 



(9) 



where we have used the classical notations of stellar 
structure theory: Hp is the pressure scale height P/pg, 
V = <91nT/91nP the logarithmic temperature gradient, 
S = — (d In p/d In T)p^ and ip = (din p/dln p,)p^- We al- 
low for a molecular weight gradient = dln/x/dlnP, 
whose contribution may be important in some cases. Let 
us also recall that, due to convective penetration, N starts 
at a finite value N c at the top of the radiation zone (Zahn 
1991). 

We introduce the vertical wavenumber k v : 



1 



'+!) 



(10) 



and observe that rk v 3> 1 for a <C N. Therefore the dif- 
ferential equation (8) may be solved by the WKB method, 
which yields the following result: 



C r 2 p 2 -^r - 1 
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(11) 



with r c designating the base of the convection zone and 
the constant C fixing the amplitude (cf. Press 1981). It 
describes a "monochromatic" wave of spherical order t 
and local frequency a, which propagates 1 with the phase 
velocity (—<r/k V} 0, r sin 9 a/m), m being the azimuthal 
order; its vertical group velocity is given by 



do- 



V ° = -dk 



a N [ 

ky 



N 2 



(12) 



Making use of the continuity equation (4), we obtain 
similar expressions for the horizontal components of the 
velocity. In particular 



rk v _3 _i / N \ 4 
£(£ +1) y g ) sin 



P£ n (cos9) 



x cos ycrt — m((j) — fit) 



k„dr 



= m ■ 
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(13) 



U r (r,6,(j),t) = U„(r)P™ (cos 9) expi[<7t — m(<p — Clt)], (6) 1 The most general solution would include a stationarywavo. 
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3. Fluxes associated with a monochromatic wave 4. Spectral distribution 



The horizontal average of the kinetic energy density is 
easily deduced from the WKB solution above: 



^p<u 2 > = -L / J 1 p{u 2 r + u 2 g + u%) sm6d6A<P 



1 n 2 i rr 

2 ^4^77 



p u r sin 6 d8 d<p 



(14) 
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Multiplying by the group velocity (12), we get the average 
flux of kinetic energy transported by a traveling wave: 



1 2 O N 2 — G 2 

Tk = 2 P<U > K^ 2 - 

. ,1 2 a 2 {N 2 -a 2 Y 
= S1 gn(M -p<u > — _ 



(15) 



with kh = ^/^(^ + l)/r being the horizontal wavenumber. 

Likewise we evaluate the mean flux of angular momen- 
tum carried by the wave: 

Tj = — f [ prsm9uj,u r sm9d9d(j) 

47T ./ ./ 



rn 



rw^hJh™ 9 * ** (16) 

(17) 



a 



where we have used (10), (13) and (15). Inserting (11) into 
(16), we verify that the angular momentum is conserved 
in this adiabatic limit 2 : 



ArKr 2 Tj{r) = Cj = est. 



(18) 



In contrast, the kinetic energy of the wave varies with 
depth, since the local frequency a depends on the rotation 
rate f2(r), but the conservation of energy is ensured in the 
inertial frame, as was explained by Bretherton (1969a) in 
the plane-parallel case. 

When radiative damping is taken into account in the 
quasi- adiabatic limit, the wave amplitude is multiplied by 
an attenuation factor exp(— r/2), where r, which is similar 
to an optical depth, is the integral 



NN 2 



3 r c ni 
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N 2 



dr 



(19) 



with K being the thermal diffusivity (see Appendix B). 
Then the angular momentum luminosity is no longer con- 
served, but decreases as 



Cj(r) = Cj(r c ) exp[-r(r)], 



(20) 



which means that some angular momentum will be ex- 
tracted from the radiative interior. 



This condition is not fulfilled in Schatzman (1993). 



From now on, we focus only on those waves which propa- 
gate towards the radiative interior and thus carry energy 
upwards (k v > 0). As we have seen above in (15), the av- 
erage kinetic energy flux transported by such a wave, at 
the top of the radiation zone (r = r c ), is given by 



LU 2 



(21) 



where u(w, kh) is the r.m.s. of the wave velocity over the 
horizontal surface, N c the buoyancy frequency, and w the 
frequency of the wave in the frame rotating with the angu- 
lar velocity f2 c of the convection zone (Press 1981; Garcia 
Lopez & Spruit 1991). 

We now turn to the evaluation of the flux carried by 
the whole spectrum of waves generated by the convective 
motions. In all likelihood the excitation of these internal 
waves takes place very close to the base of the convec- 
tion zone. There are two reasons for that. First the grav- 
ity waves do not penetrate far into the convection zone, 
since they are evanescent there with the lapse rate k h . 
And second the interface is rather sharp, with a jump in 
the buoyancy frequency and plumes penetrating into the 
radiative interior, which favors the generation of waves as 
one observes both in the laboratory (Townsend 1958) and 
in numerical simulations (Hurlburt et al. 1986, 1994). 

To obtain a crude estimate of the wave energy, we fol- 
low closely Garcia Lopez & Spruit (1991). We match the 
pressure fluctuation in the wave with that of the turbulent 
convection, thus 



u 2 (lo) = v 2 (u>), 



(22) 



where v{uu) designates the r.m.s. convective velocity at the 
frequency u>. We further take into account that convec- 
tive eddies of wavenumber k(u>) generate also waves at 
lower (horizontal) wavenumber kh- Summing the eddies 
which participate in the stochastic excitation of a wave of 
wavenumber kh, we get 



u 2 {uj,k h ) 



k{uj) 



2f, ,\ — 



= 2 1 v' 
Jo 



' dk h 



(23) 



Next we assume that the kinetic energy spectrum of 
the convective motions is represented to first approxima- 
tion by the Kolmogorov law: 
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(24) 



where v c and to c (and k c below) characterize the largest 
convective eddies. We invoke once more Kolmogorov's law 
to replace k 2 {u>) in (23) by k 2 (uj/Lu c ) 3 , and reach the fol- 
lowing expression for the flux of kinetic energy at the top 
of the radiative interior: 



Tk(t c ) = p c v\ 



■ QJc f Nc duj 



2 { l dl 
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From now on we use the spherical harmonic number £ = 
rkh, and pretend that it varies continuously from to its 
upper bound £ u (the convective scale): 



< £ < 



with £,. 



(26) 



In the limit uj c <C N c the flux integrated over the whole 
spectrum amounts to 



(27) 



a result which is similar to that given by Press for his 
monochromatic flux [cf. his eq. (90)]. Note however that 
our expression (25) for Tk vanishes when uo c — > N c , 
whereas his is singular there. 

To proceed with the evaluation of the angular momen- 
tum flux, we need to know how the spectral energy of the 
internal waves is distributed over the azimuthal order to, 
for given £. For simplicity, we shall assume that this dis- 
tribution is uniform, although it is quite possible that the 
Coriolis force, which we have neglected here, causes an 
unbalance between prograde and retrograde waves, hence 
between positive and negative to. With our simplifying 
assumption 



t r \ 3 



(28) 



d£ 



dm 
~2t' 



Referring back to (17) and (20), we obtain the following 
expression for the luminosity of angular momentum inte- 
grated over the whole wave spectrum: 



Cj(r) = Aur 



2 Pc v c 

N, 



7 

J U) c 



do,/ 



C^C 

N? 



d£ 



(29) 

exp(— r) to dm. 



In a stationary regime, one would have to subtract from 
this luminosity that, of opposite sign, which is carried by 
the waves reflected at the center, including the appropriate 
damping. 

5. The angular momentum flux 

In a non-rotating star, there is no net flux of angular mo- 
mentum, since two waves propagating in the same direc- 
tion but which are of opposite azimuthal order to have the 
same amplitude; therefore they carry (and deposit) equal 
amounts of angular momentum, which cancel each other 
because they are of opposite sign. 

The situation is different in a rotating star, because 
the Coriolis force may introduce a bias between waves of 
opposite to, and this could generate differential rotation as 



the waves undergo radiative damping. We do not consider 
this possibility here, and assume rather that differential 
rotation exists for another reason, such as mass loss or 
meridian circulation. In this case, waves of opposite az- 
imuthal order experience different damping, and the net 
effect will be an extraction of angular momentum from the 
radiation zone. When the damping is slight, the associated 
flux of angular momentum will be small, and we shall ne- 
glect it in this first approach. Instead, we concentrate on 
those waves which are completely damped because their 
local frequency vanishes at some critical depth. A similar 
damping process plays an important role in the tidal syn- 
chronization of massive binaries, as was pointed out by 
Goldrcich & Nicholson (1989). 

The local frequency may be written 



a(r) = cj — to [fl(r) — Q c ] = w — m AO(r) 



(30) 



with uj being the frequency of the monochromatic wave 
when it is emitted at the base of the convection zone, 
which rotates at the angular velocity f2 c . If the rotation 
speed Q(r) increases with depth, a will have a node at 
some location r*(u>, to) for sufficiently large (and positive) 
to. There k v — * oo and radiative damping will increase 
dramatically - formally r — > oo, but the quasi-adiabatic 
approximation leading to (19) is no longer valid. Thus at 
the critical layer r = r* this wave will deposit whatever re- 
mains from its (negative) angular momentum. The result 
is similar if the rotation rate decreases with depth, but it 
is then of opposite sign; in each case radiative dissipation 
acts to reduce the existing differential rotation. 

From now on, we shall consider only those waves which 
experience little damping before they reach the critical 
level r* where a — > 0; according to (19), their initial fre- 
quency w is such that r <, 1, or 



lu a £ I(r*)£ 3 where J(r) 



K N N] 



dr 



t r 3 ' 



(31) 



A wave which does not meet this condition is damped be- 
fore, but its contribution to the transport of angular mo- 
mentum is compensated by that of its partner of opposite 
to and we shall ignore it in this first approach. 

A final, admittedly drastic simplification will enable 
us to calculate the triple integral (29): we shall ap- 
proximate the damping exponential by the step function 
7i{r — r*). We thus take exp[— r(r)] = 1 for u > to Af2 and 
exp[— r(r)] — for u < mAil. In other words, exp[— r] 
will be removed from (29), and the integration domain in 
to will be delimited by 

< to < mm[e,(j/AQ]. (32) 
Let us recall the other limits; according to (31) and (26) 
0<£<min[(uj 4 /l)i,£ c (uj/uj c )i} (33) 
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Fig. 1. Variation with depth of the damping integral defined 
in (31), scaled by the fourth power of the convective turnover 
frequency: I/ui*. The dashed line refers to the Sun at 200 Myr, 
the continuous line to the present Sun. The radial cooordinate 
is normalized by the actual, time-dependent radius of the Sun. 



and for the frequency we have 



LO c < UJ < N c . 



(34) 



The integration can now be performed, and its results 
are reported in Appendix A. The angular momentum flux 
is an implicit function of r, through AQ(r) and I(r), whose 
variation with depth is displayed in Fig. 1. Its general ex- 
pression is rather intricate, but it simplifies in some pa- 
rameter domains. For instance, let us take the interval in 
/ which is the most relevant for the solar interior, namely 
where 



< 



< 1 



(35) 



and that in AQ(r) defined by 

/ 



UJc_ ( _Z_ 



UJc 



(36) 



(for the physical meaning of these limits we refer the 
reader to Appendix A.) In this domain the flux depends 
linearly on the differential rotation: 



PcVg 



3 (u_ 

4 \ I 



4\ 3 



1 

3 { I J to c 



(37) 



The term independent of Aft represents the contribution 
of waves which do not meet the singularity at a = 0, and 
consistent with the approximation made above, we shall 
ignore it. We thus write the result as 



Cj(r) = Cj(r c ) - 



4nr 2 p c vl 
3 NX 



An 



(38) 



6. Efficiency of the angular momentum transport 

If angular momentum is transported only by the inter- 
nal waves we have considered here, the angular velocity 
evolves with time according to 



d 

fir 2 sin 2 9 pr 2 sin 9 d6 d0 = - — C j (r) 



or equivalently 
9 ( 4o\ 1 Pc v c 9 



I 



dr 
ACL 



(39) 



(40) 



To evaluate the efficiency of this transport, we neglect 
all variations except those of fl; hence 



dfl dfl 

~8t ~ Vw ~dr~' 

whose formal solution is 

Cl(r,t) = F 



(41) 



(42) 



Within this approximation the rotation profile, whose 
slope depends on the rate of angular momentum loss 
through the wind, propagates inwards with the velocity 
V w , which is directly related to the damping integral I(r) 
defined in (31): 



V w = 



1 PcVg 1 

2 pr 2 N c uj c I c 



(43) 



Approximating p c v^ by 1/10 of the convective flux (see 
Cox & Giuli 1968), and this flux by Lq/Attt 2 , we obtain a 
crude estimate for the synchronization time t sync — r/V w : 



^sync ~ 60 



M G) R 2 



qP 
P 



Re- 



lic. 



N r . 



7j being the mean density. The variation with depth of 
t sync is depicted in Fig. 2; for this evaluation we took 
lu c = 2ttv c / Hp, £ c = 2irr c /Hp, and assumed a penetra- 
tion depth of Hp/ 10, which determines the value of A^ c 
(« 10~ 3 s _1 ). The solar models were built with the stellar 
evolution code CESAM (Morel 1996). 

The synchronization time depends somewhat on the 
age of the Sun, because the damping of the internal waves 
increases with time in the inner core, due to the steepening 
of the molecular weight gradient, but this effect is com- 
pensated by the increase of the convective frequency lo c . 
We see that the internal waves are able to extract angu- 
lar momentum from the solar interior in about 10 years, 
three orders of magnitude less than the present spin-down 
time. 
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0.2 0.4 0.6 



Fig. 2. Variation with depth of the synchronization time 
£ sy nc = r/V w (see eq. 44). As in Fig. 1, the dashed line refers 
to the Sun at 200 Myr, the continuous line to the present Sun. 

7. Conclusion 

We conclude that the internal gravity waves are quite ef- 
ficient in extracting angular momentum from the radia- 
tive interior of the Sun, a result anticipated by Schatz- 
man (1993). During the later stages of the solar spin- 
down, this process prevails over the other mechanisms 
which have been proposed so far, except perhaps mag- 
netic torquing (see Introduction). The efficiency of this 
transport has been confirmed independently by Kumar & 
Quataert (1996), but these authors considered the differ- 
ential damping between waves of opposite azimuthal or- 
der m, which redistributes angular momentum within the 
radiation zone with no net flux into the convective enve- 
lope. In contrast, we neglected here this effect altogether 
and focussed only on those waves which undergo complete 
damping at the critical level where their local frequency 
vanishes. We are implementing this transport in our rota- 
tional evolution codes, where it will compete with other 
mechanisms, such as meridian circulation, and we shall 
report the results in forthcoming papers. 

But we clearly need to move beyond the present ex- 
ploration, whose weak points ought to be recalled. First 
we have excluded the Coriolis force from the description 
of the internal waves, thus missing a possible unbalance 
between waves of opposite azimuthal order m. Second, we 
neglected the redistribution of angular momentum due to 
the differential damping mentioned above. Finally we used 
a rather crude recipe to couple the internal waves with 



the convective motions, which we assumed to obey Kol- 
mogorov's law; it could well be that the stochastic exci- 
tation at wavelengths exceeding the size of the convective 
eddies is overestimated, but the result is not too sensitive 
to the slope of the power spectrum. 

In spite of these imperfections, which have to be reme- 
died, we see only one interpretation of the short character- 
istic time we have found above. Namely that the transport 
of angular momentum through the internal waves which 
are emitted by the convection zone is responsible for the 
flat rotation profile which is observed in the Sun. 
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Appendix A: Expressions for the angular momen- 
tum flux 

The angular momentum luminosity is given by eq. (29): 

Cj(r)=4nr 2 ^-Q(r) 

with 

Let us recall that we have replaced the summations over 
integer m and £ by integrals over these quantities, which 
are assumed to vary continuously. We take m = 1 (and 
therefore I = 1) as lowest value when defining validity 
domains for various expressions. But we shall keep as 
lower bound for the integrals, because for small values of 
L they approximate better sums such that 

L I 

i^i 

We calculate first the angular momentum flux in the range 



which corresponds to case II in Fig. 3; it covers most of the 
existence domain of the integral I(r), as shown in Fig. 1. 
Then the upper limit of £, for given u>, is fixed by condition 
(33) 
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Fig. 3. Limits of the integration domain in the (£,lu) plane, in 
logarithmic coordinates. In case II, which is the most relevant 
to the solar interior, the damping line lu = J 1 / 4 ^ 3 / 4 cuts lu — lo c 
between 1 and £ c ; the corresponding integration domain is de- 
lineated in thick lines. 
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Fig. 4. Limits of the integration domain in the (m,£) plane. 

and the integration domain is split in two, according to 
whether l u is smaller or larger than lu/ AO, (see Fig. 4), or 
cquivalently whether lu is smaller or larger than cu' , 

LU C LO A c 

Subdomain a : lo c < lu < lu' (or £ u < lu/AQ) 

In the (to, £) plane the integration domain is a triangle: 



and therefore 



5 " = £' w te)' 



duj 

LU 



4\ 3 



LU 

LU C 

An 

LU C 



' dtu 

LU 

4^k■ 



(45) 



Subdomain b : lu' < lu < N c (or £ u > lu/AQ) 

The integration domain in (to, £) is now a trapezium; 
for the triangular tip we have 



r/ An d£ r 

Jo ~ k 



mdm = 



4 \avlj 



and in the rectangular part 

lj/AO * Jo 

Summing the two we get 
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It remains to integrate over lu; taking into account that 
lu c <C N c , we have 



Jui' 



LU, 



LU 



12 LU C 



and hence 
Q(r) =Q a + Qb 



3 

4 V I 



i fu4\ An 



(46) 



(47) 



which is the result quoted above in eq. (37). This expres- 
sion is valid provided lu c < lu' < N c , in other terms as 
long as 



LUc_ f J 



/Any 



' < Km < 



(48) 



So far, we have implicitly assumed that lu c < lu' < N c ; 
let us examine the other cases. 
For lu' > N r , or 



/An\ 3 lu c / I 



\ L0 C J < N c V^c 4 , 

Qb = and Q(r) does not depend on Afi anymore 

,4\ f 



(49) 



(50) 
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Appendix B: Effect of a molecular weight gradient 

A molecular weight gradient modifies somewhat the dy- 
namics of internal waves, and in particular their radiative 
damping. To take this into account, we follow the method 
adopted by Press (1981), and refer to his equations as 
(P . . ). 

The equation of motion is still (P 27) 



2 „/ 



iaV (pu r ) = gk h p 



(56) 



but the equation of state now includes the Eulerian per- 
turbation of the molecular weight p : 



(57) 



Fig. 5. The regions in the (//w*,Afl/w c ) plane where the 
angular momentum flux is given respectively by expressions 
(50), (47), (52), (53), (54) and (55). Logarithmic coordinates 
have been used again, but scales are arbitrary. 



Last we turn to the case where u>' < u r , or 



< 



( 



An 



(51) 



Then Q a = 0, and the integration of must be performed 
over the domain min(w c ,A51) < u> < N c . The result is 
(again with uo c <C N c ) 



5 1/1 AO 
In — + - In 

12 6 lu* 2 uj c 



for A£l < lo c , and 

3 



«->-(&)■ 

for Afl > u) r . 



5 1 / 
12 " 6^ 



2 , AO 
-In 

3 u c 



(52) 



(53) 



Finally we have to consider the case I of Fig. 3, where 
I > uj a , which applies to the innermost core of the evolved 
Sun (see Fig. 1). Performing the integrations as above, one 
finds 



Q(r) =Qa + Qb 



for (Afi/w c ) 4 < I/v*, and 



An 

OJr. 



Q(r) 



<4Y (^l\ 2 
i / \An) 



5 ii i , An 

In — + - In 

12 8 ujf. 2 uj c 



(54) 



(55) 



when (AO/w c ) 4 > I/lo*. 

It is easy to check that Q(r) given successively by (50), 
(47), (52), (53), (54) and (55) is a continuous function of 
An(r) and I(r), and hence of r. In Fig. 5, we display the 
domains where each of these expressions applies. 



where we have ignored the pressure perturbation, as al- 
lowed for internal waves (anelastic approximation). The 
fluctuation of p is given by 



• P V M 

lO— = -Tf-U r , 
P Hp 



(58) 



which expresses the conservation of molecular weight in 
the wave motion. 

To eliminate T', we call the heat equation 



^ = -(V ad -V)|^ + ^V 2 T', 



(59) 



K being the thermal diffusivity; using (57) and (58), and 
neglecting the vertical variation of all mean quantities 
compared to that of T", we get 



ia— = u r - K\7 2 5— 

9 9 \ T 

where 

gs 



(60) 



N-=N? + Nl = j-(V &A -V) + ^». 



H 



Next we draw (T'/T) from the equation of state (57) 
to obtain 



2N . iV 2 
(ia — K\ )p = — — pu r 



K 



ia 



Nf, 



-pu r 



(61) 



which is (P 23) plus the extra term in N 2 . Combining (56) 
and (61), we are led to 



V " + k V-^-) !>"■ ( V 2 + kl^-)pu r = 0,(62) 



.K. 
a 



which again differs from (P 28) by this term in N 2 . 

This equation is readily transformed into the disper- 
sion relation 
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- 1 



) 



(7 



. K 




which yields the following expression for the vertical 
wavenumber, in the quasi-adiabatic limit: 



and hence the value of the damping integral r(r) quoted 
above in (19): 
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